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ABSTRACT 

Pulse  mode  sonar  operation  is  analyzed  under  the 
physical  hypotheses  that 

. The  medium  is  a stationary  homogeneous  fluid. 

. Both  the  sonar  systan  and  the  scattering  objects 
are  stationary. 

. The  scattering  objects  are  rigid  bodies. 

. The  scattering  objects  lie  in  the  far  fields  of  the 
transmitter  and  receiver. 

It  is  shown  that  if  the  sonar  signal  waveform  in  the  far  field  is 


s(  X 

t,0) 

X 

then  the  sonar  echo  waveform  in  the  far  field  is 

X = (xle 

T^(a)6,a)6')s(uj,0')d6’dw  • . 

■ 

Here  s(a),0)  is  the  Fourier  transform  of  s(t,0)  with  respect  to 
T and  T^(a)0,(i)0')  is  the  differential  scattering  cross  section  of 
the  scattering  objects. 


eC|x|  - t,0) 
Ixl 


\diere 
e(T,0)  = Re^ 


^ u)  expCiro))  j 


§1.  INTRODUCTION.  This  paper  deals  with  pulse  mode  sonar  echo  prediction;  that 
is,  the  calculation  of  sonar  echoes  vdien  the  characteristics  of  the  transmitter, 
scattering  objects  and  ambient  medium  are  known.  The  physical  hypotheses  are 

• The  medium  is  a stationaiy  homogeneous  unlimited  fluid. 

• Both  the  sonar  system  and  the  scattering  objects  are  stationary. 

• The  scattering  objects  are  rigid  bodies. 

The  analysis  is  based  on  the  theory  of  scattering  for  the  wave  equation  developed 
in  the  author’s  monograph  [5].  The  principal  result  of  this  paper  is  an  asyn^stotic 
calculation  of  the  echo  waveform,  valid  when  both  transmitter  and  receiver  are  in 
the  far -field  of  the  scatterers.  The  results  show  that,  in  this  approximation, 
the  dependence  of  the  echo  waveform  on  the  scatterers  is  determined  by  the  S-matrix 
of  scattering  theory.  The  work  is  a sequel  to  the  author's  article  [6]  \diere 
similar  results  are  derived  for  plane-wave  signals. 

THH  BOUNDARY  VALUF:  PROBLEM.  In  what  follows  x = (Xj^,X2,Xj)  € denotes  the 
coordinates  of  a Cartesian  system  fixed  in  the  medium  and  t G R denotes  a time 
coordinate.  The  acoustic  field  is  characterized  by  a real-valued  acoustic  potential 
u(t,x)  which  is  a solution  of  the  wave  equation 

'y 

(i.n  -7  - v^u  = f(t,xj 

at 

The  function  f , which  characterizes  the  transmitter,  will  be  called  the  source 
function.  The  scattering  of  a single  pulse  of  duration  T , emitted  by  a 
transmitter  localized  near  a point  Xq  , will  be  analyzed.  Hence,  the  space-time 
support  of  f is  assumed  to  satisfy 

(1.2)  supp  f c {(t,x) I t^  < t < t^  + T and  |x-Xq|  < 6^} 

where  6^  and  t^^  are  constants.  The  scatterers  are  represented  by  a closed 


-2- 

f 

T T 

bounded  set  r c R with  ccmplement  n = R -r  . The  conmon  boundary  3r  = 
is  assumed  to  be  a smooth  surface.  It  will  be  convenient  to  let  the  origin  of 
coordinates  lie  in  r and 

(1.3)  r C {x|  |x|  < 6} 

It  is  also  assumed  that 

(1.4)  6 + 6o  < Ixpl 

(the  transmitter  and  scatterers  are  disjoint)  and 

(1.5)  T<|xq|-6-6q 

(the  sources  cease  acting  before  the  signal  reaches  the  scatterers). 

The  total  acoustic  field  produced  by  the  transmitter  in  the  presence  of  the 
scatterers  T is  the  solution  u(t,x)  of  (1.1)  in  Rxft  that  satisfies  the 
boundary  condition 


(1.6) 


V . 7U  - 0 
9v 


for  (t,x)  e Rx3n  , where  v is  a normal  to  , and  the  initial  condition 
(1.7)  u(t,x)  = 0 for  t < tp  and  x e n 

The  corresponding  signal  field  UQ(t,x)  , generated  by  f(t,x)  when  no  scatterers 
are  present,  is  given  by  the  retarded  potential 


(1.8) 


■ if  f , ^ 

r3 


f(t- 

x-x* 

x-x* 

dx' 


\diere  dx*  = dxjdx^dx^  . The  sonar  echo  Ug(t,x)  produced  by  the  source  function 


.1 
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f and  the  scatterers  r is  defined  by 


(1.9) 


Ug(t,x)  = u(t,x) -UgCt.x)  , tsR  , xe  n 


Both  supp  UQ(t»*)  and  supp  u(t,‘)  are  contained  in  {x|  |x-Xq|  < t-t^  + 5g^ 
vdiich  is  disjoint  from  r for  t-tg  ^ follows  that 


(1.10)  Ug(t,x)  = 0 for  t < tg  + |xg|-<S-5g  and  xe  n 

The  goal  of  this  paper  is  to  calculate  Ug(t,x)  , especially  in  the  far  field 
(|xl  » 1)  and  to  analyze  its  dependence  on  the  source  function  f and  the 
scatterers  F . 


§2.  THE  WAVE  OPERATORS  AND  PJLSE  MODE  SONAR  ECHO  STRUCTURE.  The  starting  point 
for  the  calculation  of  Ug(t,x)  below  is  a construction  of  uCt,x)  in  the  Hilbert 
space  L2(5^)  . To  describe  it  let 

(2.1)  L2(Q)  = {u(x)ln“u(x)  € for  0 < la]  < m} 

I I “l  “2  “3 

where  a = (0^,02,02)  , loi|  = aj+a2'*'U2  and  = 3'^V<iXj  8X2  3Xj  . Then  the 
operator  A : L2(n)  ■*  L2(^2)  defined  by 

(2.2)  D(A)  = 1.2(51)  {ul3u/3v  = 0 on  351}  , 

(2.3)  Au  = -V"u  for  all  u e D(A) 

is  selfadjoint  and  non-negative;  see  [5]  for  details.  The  solution  of  the  initial- 
boundary value  problem  (1.1),  (1.6),  (1.7)  is  given  by  Duhamel's  integral  [6]: 

t 

u(t,*1  = I {A  sifi(t- T)A^^^}f(T,  Odx  , t > tg 


(2.4) 


In  particular,  for  t > + T 

tfl  + T 

(2.5)  uCt,*)  = [ sin(t-T)A^^^}f(T,*)dT  = Re{v(t,*)> 

\Aere 

t +T 

(2.6)  v(t,*)  - i j ° A'^/W-i(t-T)A^/2}f(-r,-)dT  = exp{-itA^/^}h 


and 


(2.7) 


to-T 


h » i A’^^^  exp{iTA^'^^}f(T,*)dT 


In  the  special  case  where  0 = (no  scatterer)  the  operator  A vdll  be  denoted 
by  Aq  . Thus  Aq  : L2(R^)  , defined  by  D(Ag)  = and  AqU  * -Au 

is  selfadjoint  in  L23^^)  and  the  signal  UQ(t,x)  is  given  by 

(2.8)  " Re{Vg(t,0}  , t > tQ  + T 
where 

(2.9)  “ exp{-itAj^^}hQ 


and 

(2.10) 


To  conpare  u(t,x)  and  UQ(t,x)  introduce  the  operator  J : L2(n) 


defined  by 
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(2.11) 


Jh(x) 


j (x)h(x)  for  X € 


for  X e R -fi 


v^ere  j 6 c"°(R'J  , 0 < j (x)  ^ 1 , j (x)  = 1 for  |x|  > 6 and  j (x)  = 0 in  a 
neighborhood  of  F . J is  a bounded  operator  with  bound  Ijll  = 1 . It  will  be 
convenient  to  extend  the  definition  (1.9)  of  u^  by  defining  Ug(t,x)  = Re{Vg(t,x)} 
where 


(2.12)  Vg(t,x)  = Jv(t,x) -Vp(t,x)  , teR  , xeR 


The  calculation  of  the  far  field  form  of  Ug(t,x)  will  be  based  on  the 
theory  of  wave  operators  as  developed  in  [5].  The  wave  operators  and  W_  are 
defined  by  the  strong  limits 


(2.13) 


® ^ expIitA^^}.!  exp{-itA^^^} 


It  is  shown  in  [5]  that  these  limits  exist  and  define  unitary  operators 
; L-,(n)  ->•  L2(F‘^)  . It  follows  that  for  each  h e L2(fi) 

(2.14)  Jvft,-)  = J exp{-itA^'^^}h  = exp{ -itAj'^^}  W^h  + o^fl)  , t ->  + » 

•z 

where  o^(ij  denotes  an  L2(IR  ) -valued  function  of  t that  tends  to  zero  in 
L2(R^)  when  t . (2.14)  and  (2.9)  imply  that 


(2.15) 


Vg(t,‘)  = expf-itAy^}(W^h-hQ)  + o^(l)  . ^ 


This  result  is  used  below  to  calculate  the  far  field  form  of  Ug(t,x)  = Re{Vg{t,x)} 


_J 
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§3.  'IHE  FAR  FIELD  APPROXIMATIOJ  AND  THE  SCATTERING  OPERATOR.  The  scattering 

3 

operator  for  the  scatterer  r is  the  unitary  operator  S in  L2^  ) defined  by 

(3.1)  S = 

vdiere  W*  denotes  the  adjoint  of  W . A connection  between  S and  the 
approximation  (2.15)  will  be  derived  by  calculating  the  relationship  between  h 
and  hjj  . Equation  (1.10)  in^lies  that  Vg(t,x)  = 0 for  t^  + T<t<tQ  + |Xq| -6-6q 
and  X e R . It  will  be  convenient  to  choose 

(3.2)  tg  = -Ixgl  +50^*^ 

so  that  the  arrival  time  of  the  signal  at  T is  non-negative ; see  (1.10).  With 
this  convention 

(3.3)  J exp{-itA^/^}h  = exp{-itAj/^}ho  for  t^  < t < 0 
idiere 

(3.4)  tj  = tg  + T = -|Xgl  + 6g  + 6 + T 

Taking  t = 0 in  (3.3)  gives  Jh  = hg  vhile  taking  t = t^  gives 

(3.5)  expdtj^A^^^U  exp{-it^A^^^}h  = hg 

The  scatterer  T is  in  the  far  field  of  the  transmitter  if  |xq|  » 1 or,  by  (3.4), 
tj  « -1  . Combining  this  with  (3.5)  and  the  definition  (2.13)  gives 

(3.6)  hg  = W_h  + O^^(l)  , |Xgl  CO 


» 


r 

i 

*. 
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where  o (I")  is  an  L- -valued  function  of  that  tends  to  zero  in 

Xq  ^ u z 

when  |Xq|  ^ 00  . Multiplying  (3.6)  by  S gives 

(3.7)  W^h  = Shp  + o^^(l)  , IxqI  - oo 

because  W*W  = 1 (W  is  unitary).  Combining  (3.7)  and  (2.15)  gives 

(3.8)  ^ exp{-itAy^}(S-l)hQ  o^(l)+Ojj^(l) 

Note  that  the  term  o (1)  tends  to  zero  in  L-(R  ) when  |x-|  -+-00  uniformly 

Xq  z u 

in  t because  exp{ -itAj'^^}  is  unitary.  Equation  (3.8)  shows  that,  in  the  far 
field  approximation,  the  dependence  of  the  echo  waveform  on  the  scatterer  is 
determined  by  the  scattering  operator. 

The  approximation  (3.8)  is  used  in  56  to  derive  an  explicit  integral  formula 
for  the  far  field  echo  waveform.  The  derivation  is  based  on  a known  integral 
representation  for  S , formulated  in  §4,  and  the  theory  of  asymptotic  wave 
functions  of  [5]  which  is  sunmarized  in  §5. 

§4.  THE  STRUCTURE  OF  THE  SCATTERING  OPERATOR.  The  steady- state  theory  of 
scattering  and  associated  eigenfunction  expansions  for  A are  reviewed  briefly 
in  this  section  and  applied  to  the  construction  of  the  scattering  operator  for 

r . 

Aq  is  a selfadjoint  operator  in  ® purely  continuous  spectrum 

and  the  plane  waves 


f 


(4.1)  Wq(x,p)  = (Zir)  exp{ix*p}  , p e 


are  a complete  family  of  generalized  eigenfunctions.  The  corresponding  eigen- 
; function  expansion  is  the  well-known  Plancherel  theory  of  the  Fourier  transform 

i (see  [5],  Ch.  6). 

i 

I 

Generalizations  of  the  Plancherel  theory  to  acoustic  scattering  by  bounded 
objects  were  first  given  by  N.  A.  Shenk  [2]  and  Y.  Shizuta  [4].  In  this  work 
the  generalized  eigenfunctions  are  the  distorted  plane  waves 


M 


1 


(4.2) 


s 3 

w+Cx,p)  = w-(x,p)  + w^(x,p)  , X e n , p e R 


which  are  characterized  by  the  properties  that  w^(x,p)  is  locally  in  D(A) 
(i.e.,  (J)w^(*,p)  e D(A)  for  all  (|)  e c“(R^))  , 

(4.3)  (V^+|p|^)w^(x,p)  = 0 for  xen 


(4.4) 


3j7r'i|p|wt-o(-j^) , ixi 


For  the  existence,  uniqueness  and  construction  of  w^(x,p)  see  [2, 3, 4, 5, 6]. 
Physically,  w^(x,p)  is  the  steady-state  scattered  field  produced  when  the 
plane  wave  (4.1)  is  scattered  by  F . It  has  the  far  field  form  [5,6] 


(4.5) 


±i|p||x|  -1 

^ 4'tr|x| — ^ 


where  0 = x/|x|  . T^(p,p')  , the  scattering  amplitude  or  differential  scattering 

cross  section  of  F , is  defined  for  all  p and  p'  in  R^  such  that  jpj  * lp'| 
and  has  the  symmetry  properties 
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4 

t 

\ ' 

h 


i ! 


(4.6) 


T+CP.P')  = TJ-P',-P)  = T^rp,-W  = T^C-P.'P*) 


vdiere  the  bar  denotes  the  complex  conjugate. 

The  connection  between  S and  T^(p,p')  is  based  on  the  eigenfunction 
expansion  theorem  for  A . The  latter  states  that  for  all  h e L2(S^)  the 
limits 


(4.7)  h^(p)  = (4^h)(p)  = L^CR  )-lim 

exist,  where  = n {x|  jxj  < M)  , and 


w^(x,p)l  h(x)dx 


(4.8) 


h(x)  = L2(J^)-lim  j w^(x,p)h^(p)dp 
M + 00  I P 1 M 


Moreover,  the  operators  ; 1^2 (^2)  L2(R^)  are  unitary  and  for  each  ^xjunded 

measurable  function  ’t'(A)  on  \ > 0 


(4.9) 


(4>^4'(A)h)(p)  = M'(|pl“)4>^h(p) 


The  Fourier  transform  will  be  denoted  by 


(4.10) 


h(p)  = ($h)(p)  = L2(R^)-lim 

M-»-<»  |x|  <M 


Wq(x,^  h(x)dx 


An  exposition  of  the  eigenfunction  expansion  theory  is  given  in  [5].  In  what 
follows  the  relations  (4.7),  (4,8)  are  written 


(4.11) 


h+(p)  = (x,p)  h(x)dx 


Q 
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C4.12) 


hCx)  = 


w^(x,p)h^(p)dp 


IT 


for  brevity.  However,  the  integrals  are  not  convergent,  in  general, and  (4.11), 

(4.12)  must  be  interpreted  in  the  sense  of  (4.7),  (4.8). 


The  wave  operators  defined  by  (2.3)  are  known  to  have  the  representations 
[5] 


(4.13) 


= 4'*4>_  , W_  = 4*4^ 


Combining  (4.13)  and  (3.1)  gives  the  representation 

(4.14)  S = W^W*  = 4*5'}' 


vdiere  the  operator 

(4.15)  S = 

/s 

is  called  the  S-matrix  for  the  scatterer  r . The  operator  S 1 has  the  inte- 
gral representation 


(4.16) 


(S-l)h(p) 


■M 


2(2tt) 


17? 


T^(P,|p|0')h(|p|0')de' 


vdiose  kernel  the  differential  scattering  cross  section  of  F . The  integration 

2 3 

in  (4.16)  is  over  the  points  0*  of  the  unit  sphere  S in  F . The  first 
proof  of  (4.16)  for  acoustic  scattering  is  due  to  Shenk  [2]. 


§5.  MODE  SONAR  SIGNALS  IN  TtiE  FAR  FIELD.  The  signals  UQ(t,x)  origi- 

nate in  the  region  Ix-x^j  < 6q  and  reach  points  x in  the  far  field,  charac- 


* 
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terized  by  |x-Xq|  » 1 , after  a time  interval  of  magnitude  conparable  with 
|x-Xq|  . Hence  the  far  field  form  of  u^Ct.x)  coincides  with  its  asymptotic 
form  for  large  t . The  latter  is  provided  by  the  theory  of  asynptotic  wave 
functions  developed  in  [5].  The  theory  is  applied  here  to  determine  the  far 
field  form  of  UpCt.x)  . 

The  complex  wave  function  v^ft.x)  defined  by  (2.9),  (2.10)  has  the  Fourier 


representation 


(5.1) 


Vn(t,x)  = (2ti) 


■5/2  ■ 


exp{i(x-p-t|p|)}h„(p)dp 


Equations  (2.10)  and  (4.9)  imply  that 


(5.2) 


where 


hnfP^  = (2Tr)^^^  ilpl’^  f(-|pl,p) 


(5.3)  ?(w,p)  = (2Tr)  ^ exp{-i(ta)+x*p)  }f(t,x)dtdx 


is  the  4-dimensional  Fourier  transform  of  f . Note  that  (5.2)  suggests  the 
concept  of  a non-radiating  source  function,  f is  said  to  be  non-radiating  if 


(5.4) 


2 2 2 A 

f(t,x)  = 3 u„/3t^-v'^u„  , supp  u-  bounded  in  k 


In  this  case  Up,(oj,p)  exists  and  (5.4)  is  equivalent  to 


(5.5) 


f(w,p)  = (|p|^-(ij^)uQ(u,p) 


Equations  (5.2)  and  (5.5)  imply  fif.(p)  = 0 and  hence  u»(t,x)  = 0 for 


t > to  ^ T . 
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The  asyn^Jtotic  wave  function  associated  with  u^Ct.x)  = ReiVpCt.x)}  is 


defined  by 


(5.6) 


u~(t,x)  = s(|xl-t,6)/|x|  , X = |x|e 


where  s e L2CR  x S^)  is  defined  by  (see  [5] , Ch.  2) 


(5.7) 


00 

i(T,0)  = Re  |(2tt)'^^^  I exp{iTw}(-ia))i^((De)du)| 

0 

00 

= Re  I exp{iT(ij}f (-aj,a)0)da)| 


Direct  calculation  of  s(t,0)  fran  (1.8)  yields  the  alternative  representation 


(5.8) 


s(t,0)  = ^ f(-T+0*x,x)dx 


It  was  shown  in  [5]  that  u"  describes  the  asynptotic  behavior  of  Uq  in 
L2(JI^)  for  t -»■<»: 


(5.9) 


UQ(t,0  = u“(t,*)  + 0^(1)  , t -► 


The  integral 


(5.10)  E(u,K,t)  =jj  {lVu(t,x)|^  + (3u(t,x)/3t)^}dx 

K 

may  be  interpreted  as  the  acoustic  energy  in  the  set  K c at  time  t . 

It  was  shown  in  [5]  that  if  h^  e L^CR^)  then  u“(t,x)  converges  to  UQ(t,x) 
in  energy  when  t-*^"  . More  precisely, 
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(5.11)  3uQ(t,x)/3x^  = UpjCt.x)  +0^(1)  , t -►  <»  , j = 0,1, 2, 3 

v^iere  Xq  = t , 


(5.12) 


UQ,j(t,x)  = Sj(lx|-t,0)/|xl 


(5.13) 

and 


Sq(t,0)  = -3s(t,0)/9t 


(5.14) 


Sj(T,0)  = -6jSQ(T,0)  for  j = 1,2,3 


The  result  (5.11)  was  used  in  [5]  to  calculate  the  asymptotic  distribution 
of  energy  in  cones 


(5.15) 


C = {x  = r0|  r > 0 , 0 e Cq  c 


Applying  the  results  to  the  signal  UQ(t,x)  generated  by  f(t,x)  gives 


(5.16) 


E(Uj^,C,<»)  = lim  r:(Un,C,t)  = it 


0’ 


lf(-lpl .P)l  dp 


In  particular,  the  total  signal  energy  introduced  by  the  source  function  f is 


(5.17) 


E(UgJR  ,«)  = TT 


|f(-|p|.P)rdp 


§6.  PULSE  MODE  SONAR  ECHOES  IN  THE  FAR  FIELD.  Equation  (3.8)  implies  that 
the  echo  Ug(t,x)  satisfies 
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(0.1)  = Re{exp(-itAj/2)(S-l)hQ}  +o^(l) 

The  first  term  on  the  right  has  the  same  form  as  the  signal  but  with  h^  re- 
placed by  (S-l)hQ  . It  follows  from  the  results  of  §5  that 


(6.2) 


where 


Ug(t,x)  = u“(t,x)  +0^(1) 


(6.3) 


u"(t,x)  = e(|x|-t,e)/|x|  , x=  |x|e 


(6.4) 


5(t,0)  = Re{(2TT)~^^^  I exp(iTu))  (-i(u)  ((S-l)hQ)  ((D0)da)} 


Now  by  (4.14) 


(6.5) 


((S-Dhg)  = $(S-l)$*4hp  = (S-l)hQ 


and  hence  by  (4.16) 


• r /s 

(6.6)  ((S-l)hQ)"(u)0)  = T^(a)0,a)0')hQ(a)0')d0' 


Combining  (6.4)  and  (6.6)  gives 


GD 

(6.7)  e(T,0)  = ^ Re  I I exp(iTti))(i)^  | T_|_((D0,a)0')hQ((D0')d0'da)| 


Finally,  by  (5.2) 


(6.8)  e(T,0)  = T77  ^ I j T^(i»)0,u)0')f(-u),u)0')d6’d<jo| 
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Equation  (6.2)  implies  that  u^,  gives  the  far  field  form  of  . For 
|x»l  » 1, r is  in  the  far  field  of  the  transmitter  and  the  term  o (1)  is 

u Xq 

small,  uniformly  for  all  t . For  receivers  in  the  far  field  region  lx|»l 
for  r , the  echo  u^(t,x)  arrives  at  times  t » 1 and  hence  the  term  o^(l) 
is  small. 


§7.  CONCLUDING  RB1ARKS.  Actual  sonar  transmitters  do  not,  of  course,  generate 
signals  by  means  of  a source  function  f(t,x)  . However,  the  purpose  of  a 
well  designed  transmitter  is  to  generate  a signal  with  a prescribed  waveform 
s(t,0)  . Now 


(7.1) 


s(w,9) 


00 

exp(-ia)T)s(T  ,0)dT 

-oo 


1/2. 


f (-m,a)0) 


(see  (5.7))  and  hence 

00 

(7.2)  e(T,9)  = Re  I I I exp(iTw) 

0 


ID 


T^((De,tDe')s(u),0')d6'da) 


In  particular,  the  transmitter  characteristics  influence  the  echo  waveform  only 
through  sfT,9)  . Hence,  (7.2)  is  applicable  to  real  transmitters  with  known 
waveforms  s(t,0)  . 

It  is  known  that  T^(a)0,u)0')  is  a meromorphic  function  of  u with  poles 
in  the  lower  half  plane  [3].  The  other  functions  in  the  integrand  of  (7.2)  are 
entire  holomorphic  functions.  Hence,  the  integral  in  (7.2)  can  be  transformed 
by  deforming  the  contour  of  integration  in  the  complex  w-plane.  This  leads 
to  an  expansion  of  the  echo  waveform  of  the  type  occurring  in  the  singularity 
expansion  method  [1]. 
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IS.  AB>^RACT 

Pulse  mode  sonar  operation  is  analyzed  under  the  physical  hypotheses  that 

(the  medium  is  a stationary  homogeneous  fluid. 

Both  the  sonar  system  and  the  scattering  objects  are  stationary. 
The  scattering  objects  are  rigid  bodies,  '' 

'fhe  scattering  objects  lie  in  the  far  fields  of  the  transmitter 
and  receiv'er.  ^ 

It  is  shown  that  if  the  sonar  signal  waveform  in  the  far  field  is 


— • X = |x|0 
then  the  sonar  echo  waveform  in  the  far  field  is 

, |,|e 

where 

* 

e(r,e')  = Re' -1  ( w explirw)  T (u)0,wyjs(u),O' )d0'du)  • . 

^ JQ 

Here  s(a),0)  is  the  Fourier  transform  of  s(T,a)  with  respect  to  t and  T^((u0,u)0') 
is  the  differential  scattering  cross  section  of  the  scattering  objects. 


where 
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